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ABSTRACT 

Using the Euchdean value of < V/V^^.^ > as a cosmo logical distance indicator, 
we derive the isotropic-equivalent characteristic peak luminosity of gamma-ray 
bursts both longer and shorter than 2 s. The short bursts have essentially the 
same characteristic peak luminosity of 0.6 x 10^^ erg (0.064s)^^ as do the long 
bursts. This may apply also to bursts with durations less than 0.25 s. The local 
space density of short bursts is around three times lower than that of long bursts. 

Subject headings: cosmology: observations — gamma rays: bursts 

1. Introduction 

Since 1997, our understanding of gamma-ray bursts (GRB) has increased enormously 
through the optical identification of afterglows, determination of redshifts from optical spec- 
tra and generally from afterglow studies covering a large range of energies, from radio to 
X-rays. These studies suggest that GRBs are associated with massive stars and hence that 
the burst rate may be linked to the rate of star formation in galaxies. The number of GRBs 
so observed is still rather hmited, and those successfully studied are all long bursts. We 
use the classification of GRBs in long bursts and short bursts based on the distribution of 
durations Tgo in the BATSE data (Kouveliotou et al. 1993), which show a minimum around 
2 s. 

No afterglows have been observed so far for short bursts (Gandolfi et al. 2000). As a 
consequence, there is no direct knowledge of their redshifts, luminosities and space densities. 
There has been a suspicion that short GRBs would be at smaller redshifts and be of lower 
luminosities than long bursts, perhaps due to a remark in the abstract of a paper by Tavani 
(1998) that short bursts show little if any deviation from a Euclidean distribution but this 
is not supported by data discussed in the paper. 



- 2 - 



In this Letter, we briefly discuss efforts to derive tlie luminosity function of long bursts 
and then present a derivation of the characteristic peak luminosity L* and the local space 
density for both long and short bursts. Even though the total gamma-ray energy may be 
physically of greater interest (Frail ct al. 2001), we concentrate on peak luminosities since 
the detection of GRBs is based on count rates, not on the time-integrated flux or fluence of 
the burst. 

Before any redshifts of GRBs were known, Mao, Narayan, & Piran (1994) used the 
Euclidean value of < K/V^ax > as a relative distance indicator to show that long and short 
bursts had the same peak luminosity to within a factor of 2, assuming they were standard 
candles. The redshifts that are now available for long bursts have led to the development of 
several luminosity indicators, such as the spectral lag derived from cross-correlation of two 
spectral channels (Norris, Marani & Bonnell 2000) and the variability in the time profile 
(Fenimore & Ramirez- Ruiz 2000). These luminosity indicators can in principle be used to 
derive the luminosity function of GRBs, including its evolution with redshift. Norris, Scargle 
& Bonnell (2001) have shown that the spectral lags for short bursts are much smaller than 
those of long bursts and conclude that the lag magnitude is discontinuous accross the 2 s 
valley between long and short bursts. 

We have used the Euclidean value of < V/Vra&x > as a cosmological distance indicator 
to first derive characteristic luminosities (Schmidt 1999b) and then the luminosity function 
of GRBs (Schmidt 2001). This method, which makes use of a large sample of GRBs, cir- 
cumvents the current weakness of methods based on the small number of redshifts observed 
so far. The price we pay for using < V"/V"max > as a distance indicator is that an assumption 
has to be made about the evolution of the luminosity function with redshift. We assume 
that this evolution tracks the star formation rate (SFR) on the expectation that GRBs are 
associated with massive stars. We employ the parametrizations of Porciani & Madau (2001), 
in particular their SFR model SF2 in which the co-moving space density rises by an order 
of magnitude near redshift z — 1 and then remains roughly constant for 2; > 2 (Steidel et al. 
1999). 

Our application of < V/Kiax > as a distance indicator used the BD2 sample (see Sec. 
2) which is based on BATSE DISCLA data on a time resolution of 1024 ms. The resulting 
characteristic luminosities (Schmidt 1999b) and luminosity functions (Schmidt 2001) there- 
fore applied to GRBs with a duration exceeding 1 or 2 s. In this Letter we employ the 
same method to derive characteristic luminosities L* using the BATSE catalog. By using 
the GRBs detected with a time resolution of 64 ms, we can derive the L* of short bursts. 
We will actually employ time resolutions of both 64 and 1024 ms, and obtain L* for both 
short and long bursts, so that we can directly compare them. 
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In Sec. 2, we review the application of < V/Vuiax > in deriving the luminosity function 
of long GRBs from the BD2 sample without using redshifts and recall the main results 
obtained. Following a discussion of data from the BATSE catalog in Sec. 3, we derive 
characteristic luminosities for both long and short GRBs in Sec. 4. The results are discussed 
in Sec. 5. Throughout this Letter, we use a flat cosmological model with Hq = 65 km 
Mpc-\ Qm = 0.3, and = 0.7. 

2. Luminosity Function Of Long Bursts Derived From The BD2 Sample 

The BD2 sample of GRBs is based on BATSE DISCLA data consisting of the continuous 
data stream from the eight BATSE LAD detectors in four energy channels on a timescale of 
1024 ms (Fishman et al. 1989). The sample was derived using a software trigger algorithm 
requiring an excess of at least 5cr over background in at least two detectors in the energy 
range 50—300 keV. The initial version was described in Schmidt (1999a) and a revision in Sec. 
2 of Schmidt (1999b). The BD2 sample covers a period of 5.9 yr from TJD 8365 - 10528. 
It contains 1391 GRBs of which 1013 are hsted in the BATSE catalog. The value of < 
V/Vr^^ >= 0.336 ± 0.008. The sample of 1391 GRBs effectively represents 2.003 yr of full 
sky coverage, corresponding to an annual rate of 694 GRBs. 

The derivation of the luminosity function of GRBs (Schmidt 2001) was based on the 
correlation of the Euclidean value of < V/Vmax > with spectral hardness. Given that the 
Euclidean value of < V/Vuia^ > for a well defined sample of cosmological objects is a dis- 
tance indicator, we interpreted the correlation of < V"/Vmax > with spectral hardness as a 
luminosity-hardness correlation. The luminosity function was derived as the sum of the lu- 
minosity functions of four spectral hardness classes. It can be characterized approximately as 
consisting of two power laws of slopes —0.6 and —2, respectively, with an isotropic-equivalent 
break peak luminosity of logL* ~ 51.5. The luminosity function ranges approximately from 
logL*-1.5to logL* + 1.0. 

In the derivation of the characteristic luminosity L* of GRBs from BATSE data in Sec. 
4, we will assume that the shape and extent of the luminosity function is that of the broken 
power law just described and derive the value of logL* from < V/Knax >• This is essentially 
the method used in Schmidt (1999a), where we varied the assumed extent and shape of the 
luminosity function to study the effect on the derived value of logL*. 
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3. Using Data From The BATSE Catalog 

3.1. Evaluating l^/Kiax 

We will be using < V/Vmax > values derived from the BATSE 4B catalog^ for bursts 
both longer and shorter than 2 s. Before we apply < V/Vmax > as a cosmological distance 
indicator, we consider how the individual V/Vmax values in the BATSE catalog are derived, 
and also how they are affected by the imposition of duration limits. 

The BATSE catalog lists for individual bursts the count rate in the second brightest 
illuminated detector, Cmax, as well as the minimum detectable rate C^j„. The value of 
V/Vraax is then simply derived as (Cmax/Cmin)'^^'^ ■ In contrast, the values of V/Vmax in 
the BD2 sample have been derived through simulations, in which the burst is moved out 
in Euclidean space in small steps with a corresponding reduction in its amplitude. At each 
step, the full detection algorithm is re-employed to set the background and detect the burst. 
Once the burst is not detected any more, the value of V^/Knax is simply derived from the 
reduction factor. During this process of removal, the burst may get detected later and later 
depending on the time profile. The background time window, which precedes the detection 
by a fixed time interval, may start to include some burst signal. The final detection is usually 
made on the peak of the burst but in some cases where the burst signal preceding the peak is 
high and enters the background, the final detection may be off the peak. Given that in these 
cases upon removal the burst drops out earlier than expected from Cmax/Cmin, the actual 
vahie of K/Vmax "will be larger. In practice, the effect depends much on the time profile. The 
net effect for a sample of bursts is that {Cmax/Cmin)~^^^ is an underestimate of < V/V^rnax >■ 

Next we consider the effect of imposing duration limits, such as Tgo < 2 s. In the 
simulations carried out on the BD2 sample described above, we found that the duration of 
the burst decreased as it was moved out until at the limit of detection it was 1 or 2 s. A 
qualitatively similar effect has been described as a ffuence duration bias for GRBs in the 
BATSE catalog (Hakkila et al. 2000). Suppose our sample is set by a restriction involving a 
limiting duration Tu^- In deriving V/Vmax, we should strictly apply two simultaneous limits, 
namely C^j„ and T^^, as was done in the first V/Vraax application (Schmidt 1968). However, 
evaluating the effect of T^^ on < V/V^ax > would require simulations of the derivation of 
Tgo for BATSE GRBs which are not available. Therefore, we limit ourselves to finding the 
sign of the systematic error in < V/Vuiax > if we ignore Tum- 

Consider the case of a GRB with a duration Tgo > Ti. As we move the burst out, Tgo 



See http://gamniaray.nisfc.nasa.gov/batse/grb/catalog/4b catalog maintained by W. S. Paciesas et al. 
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will decrease and may become smaller than 71 before it becomes undetectable. Therefore 
the reduction factor is smaller and V/V^ax is larger. Ignoring the lower limit Ti leads to an 
underestimate of V^/Knax- In the case of an upper limit T2, a GRB with Tqq > T2 which does 
not belong to the sample, may become part of it when its T90 becomes shorter than T2 upon 
removal. In this case, ignoring the upper limit T2 produces a V^/Knax that is an overestimate. 
It should be emphasized that these considerations are purely qualitative; the actual effects 
depend on such factors as the time profile of the burst, the way the background is set, etc. 



3.2. Effective Coverage 

In order to derive the rate of GRBs per unit volume, we need to have an estimate of 
the effective full sky coverage of the GRB sample used. For the BD2 we have evaluated the 
efficiency at 33.8% leading to an effective full sky coverage of 2.003 yr (Schmidt 1999a). The 
total sample of 1391 GRBs then corresponds to a rate of 694 GRB yr"^. 

The 4B catalog gives an annual rate of 666 GRB yr^^, presumably based on the most 
sensitive detections, which arc at the 1024 ms timescale. Considering that the S/N limit 
of the 4B catalog is 5.5a and that of the BD2 sample 5a, these rates arc quite consistent. 
We only consider GRBs detected while the BATSE on-board trigger was set for 5.5(7 over 
the energy range 50 — 300 keV. Among those, the 4B catalog contains 612 GRBs for which 
Cmax/Cmin > 1 at the 1024 ms timescale (those < 1 were detected at time scales of 64 ms 
or 256 ms). This corresponds to an effective full sky coverage for the purpose of this work 
of 0.92 yr. We assume that this value also appfies to the 64 ms timescale. 



4. Characteristic Luminosities For Long and Short Bursts 

The method used to derive the isotropic-equivalent characteristic peak luminosity L* 
for a given value of < V/Vmax > is essentially the same as that used before (Schmidt 1999b). 
Based on our discussion of long bursts in Sec. 2, the local luminosity function of peak GRB 
luminosities L, defined as the co-moving space density of GRBs in the interval logL to 
log L + d log L, is 

$o(L) = 0, for logL < logL* - 1.5, (la) 
$„(L) = c„(L/L*)-°^ for logL* - 1.5 < logL < logL*, (lb) 
$o(L) = Co{L/L*)-^-^, for logL* < logL < logL* + 1.0, (Ic) 
$o(L) = 0, for logL > logL* + 1.0. (Id) 
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We assume that the GRB rate as a function of redshift follows the SFR model SF2 (see 
Sec. 1). The median value of the spectral photon index for the (long) bursts in the BD2 
sample is —1.6 (Schmidt 2001). For the short bursts wc adopt an index of —1.1 to reflect 
their larger average hardness ratio (Kouveliotou et al. 1993). 

The median 4B limiting photon flux for 1024 ms detection in the 50 — 300 keV range, 
including the effect of atmospheric scattering, is 0.25 ph cm~^ s~^ (Mccgan et al. 1998). We 
adopt the same value for the BD2 sample for which atmospheric scattering has not yet been 
evaluated. For GRBs detected on the 64 ms timescale, the derivation of L* from < l^/Knax > 
is carried out on this timescale, so the resulting luminosity L* is produced in ergs (0.064s)~^. 
In order to convert to this system from a timescale of 1024 ms, we compared for GRBs with 
Tgo > 2 s their peak fluxes at timescales 1024 ms and 64 ms given in the BATSE catalog. 
The average 1024/64 flux ratio is 0.68, reflecting the variabihty of long GRBs at their peak 
at subsecond timescales. Accordingly, we convert peak luminosities per 1024 ms into ones 
per 64 ms by dividing by 0.68 x 16 = 10.9. 

In order to provide a comparison with past work, we first derive L* for long bursts. 
The top line of Table 1 shows the results for the BD2 sample of GRBs. The value of L* 
is 70% larger than that found from the more sophisticated derivation in Schmidt (2001), 
which involved splitting the sample in groups of different spectral hardness. The resulting 
luminosity function was not exactly a broken power law, which causes the above difference. 

The remaining entries in Table 1 are all based on data from the 4B catalog. The first 
three can be directly compared with the BD2 results since they all concern long bursts. The 
logL* values show a range of 0.3, with a systematic offset of around +0.2 from the BD2 
value. The local densities are systematically lower than that from the BD2, partly due to 
the higher luminosity. 

The next row of Table 1 gives the results for short bursts with T90 < 2 s. The charac- 
teristic peak luminosity L* is in the middle of the range given for long BATSE bursts. The 
local space density of short bursts is ~ 3 times smaller than that of long bursts. In order to 
investigate whether there might be a trend among the short bursts, we show further results 
for T90 < 0.5 s and < 0.25 s. The formal mean errors in logL* for the short bursts are ±0.20, 
±0.25, and ±0.35, respectively. We see no trend of logL* with duration among the short 
bursts. 
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5. Discussion 

The systematic effect of using Cmax/ Cmin in tlie derivation of < K/Knax > discussed 
in Sec. 3 applies to all BATSE entries in Table 1. Therefore all values of < V/Vmax > are 
systematically too small, and the resulting values of log L* too large. The actual effect is 
hard to estimate because it depends on the burst profiles and the methodology of setting the 
background. For the BD2 sample we find that using Cmax/ Cmin results in an underestimate 
of < V/Vraa^ > of 0.009 leading to an overestimate of logL* by 0.08. Considering this offset 
the agreement between the first two rows of Table 1 is satisfactory considering that the 
formal errors (derived from the mean errors of < V/Kaax >) in logL* are ±0.07 and ±0.10, 
respectively. For the BATSE bursts with Tgg > 2 s in the next two rows of Table 1, the 
systematic effect caused by the duration limit may result in an overestimate of logL*, but 
no obvious effect is evident. 

The short bursts with Tgo < 2s have a value of logL* that is entirely consistent with 
that of the long bursts. The systematic effect of the upper limit for the duration will be an 
underestimate of logL*. This effect may well be small for short bursts, but only simulations 
can tell. 

The peak luminosities and local space densities given in Table 1 are isotropic-equivalent 
values and are all based on SFR model SF2 (Sec. 1). If, instead, we use SFR models SFl or 
SF3 the derived values of logL* typically change by —0.17 and 0.14, respectively, essentially 
the same for short and long bursts. 

Assuming that the BATSE values of logL* are offset by ±0.08, we find that for our 
given assumptions about the shape and the evolution of the luminosity function, long bursts 
with Tgo > 2 s and short bursts with Tgg < 2 s have the same value of the characteristic peak 
luminosity L* ~ 0.6 x 10^^ erg (0.064s) As a consequence, short bursts will have lower 
radiated energies than long bursts, which may have a major effect on the afterglows of short 
GRBs (Panaitescu, Kumar, & Narayan 2001). Among the short bursts there is no evidence 
for any substantial change in L* for durations as short as 0.25 s. 
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Table 1. Characteristic Luminosity and Space Density From Various GRB Samplesf 



Sample 


At 


Tgo 


n 


< V/Vn,^ > 


a 


f b 
J lim 


logL*^ 


o 

Ho 


BD2 


1024 ms 




1391 


0.336±0.008 


-1.6 


0.25 


50.67 


0.339 


BATSE 


1024 ms 




612 


0.312±0.012 


-1.6 


0.25 


50.88 


0.240 


BATSE 


1024 ms 


> 2s 


469 


0.291±0.013 


-1.6 


0.25 


51.05 


0.147 


BATSE 


64 ms 


> 2s 


323 


0.376±0.016 


-1.6 


1.00 


50.72 


0.250 


BATSE 


64 ms 


< 2s 


141 


0.354±0.023 


-1.1 


1.00 


50.92 


0.075 


BATSE 


64 ms 


< 0.5s 


91 


0.361±0.029 


-1.1 


1.00 


50.86 


0.040 


BATSE 


64 ms 


< 0.25s 


49 


0.347±0.042 


-1.1 


1.00 


50.98 


0.023 



^We use Hq = 65 km s'^ Mpc-^ Vtu = 0.3, and J^a = 0.7. 

^fiim is the hmiting flux in ph cm~^ s~^ over the energy range 50 — 300 keV. 

^L* is the isotropic-equivalent characteristic peak luminosity in erg (0.064s) ~^ in 
the 50 - 300 keV band. 

^Po is the local {z — 0) GRB rate, in units of Gpc~^ yr~^- 



